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Abstract — Using the Green’s dyad technique based on cuboidal meshing, we compute the elec¬ 
tromagnetic field scattered by metal nanorods with high aspect ratio. We investigate the effect 
of the meshing shape on the numerical simulations. We observe that discretizing the object with 
cells with aspect ratios similar to the object’s aspect ratio improves the computations, without 
degrading the convergency. We also compare our numerical simulations to finite element method 
and discuss further possible improvements. 

Keywords: Dyadic Green Method; Finite element method; Metallic nanoscatterer; Surface plas- 
mon resonance; 


1. INTRODUCTION 

Plasmonic nanostructures concentrate light on subwavelength scales mm, opening the way to 
promising applications such as nano-optical antennas mm , or surface enhanced spectroscopies mm- 
In the last decade, important efforts have been done in modeling these structures [7]. Particularly, 
the electric field rapidly decays away from the metal surface so that dedicated methods have to be 
developed to correctly describe near-field optics properties of plasmonic systems. Additionnally, 
high aspect ratio shapes are often needed to tune resonance frequencies mm- Fast field decay and 
high aspect ratio lead to serious numerical difficulties for methods based on the discretization of 
the object volume, notably in term of memory cost and computational time. 

We are specifically interested in the Green’s dyad technique (GDT) (also called volume integral 
method) [10] since it naturally satisfies boundary conditions both at the object surfaces and in 
the far-field. Moreover, it easily includes a substrate supporting nanoparticles without the need 
to discretize it m- In addition, as far as the Green’s dyadic is numerically computed in the 
source (nanostructures) region, ah the electromagnetic properties of the system are easily obtained 
with practically almost no additionnal computing cost m- This concerns obviously the electric 
and magnetic fields in the very near-field of the object m but also scattering properties in the 
radiation zone nans]. Last, the local density of optical states (LDOS), an intrinsic property of 
the nanostructure, is also easily deduced from this formalism m- 

Recent works were devoted to improve the GDT when applied to metallic nanostructures. The 
GDT relies on the discretization of the scatterers. Because of strongly varying fields, the key point is 
to correctly describe the fast variation of the Green’s tensor over an elementary cell. Let us mention 
the regularization scheme proposed by Kottmann and Martin developed for 2D-elements but also 
transposable to 3D-nanostructures nzj- Instead of regularizing the Green’s tensor, Chaumet and 
coworkers directly computed its integral over the volume of cubic meshes m- To this aim, they 
considered a Weyl expansion of the tensor and performed a numerical integration in the reciprocical 
k-space. In the present work, we first transform the volume integral to a surface integral over a 
cuboidal mesh m before a numerical evaluation in direct space. This avoids the singularity at the 
mesh center and convergence difficulties. Moreover, this permits to consider more complex meshing. 
Note finally that our method differs from surface integral method (also called boundary element 
method) where the object surface is discretized, including the substrate if necessary [20l EU [22] . 
Surface integral method relies on the introduction of equivalent surface currents that depends on 
the excitation field. Differently, the volume integral method relies on a volume discretization of 
the object, excluding the substrate, but in the present work involves integration over the surface 
of the meshes. This leads to the evaluation of the Green’s dyad associated to the whole structure 
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that is independent on the illumination conditions so that various problem can be treated by post¬ 
processing. 

In the following, we define a benchmark configuration that consists of a silver nanorod. In order 
to assess the reliability of our results, we compare our numerical simulations to those obtained 
using a commercial software (COMSOL Multiphysics) based on the finite element method (FEM) 
[23] . The objective of this comparison is twofold. FEM efficiently considers complex shape but 
necessitates a careful design and positioning of a perfectly matched layer (PML) to avoid artefact 
reflexion at the boundaries of the computational window. As we will see, it becomes a very sensi¬ 
tive parameter in presence of a substrate. On the opposite, GDT easily considers nanostructures 
deposited on a substrate but is generally limited to simpler shapes. By comparing the two methods, 
we are able to estimate the error. Reciprocally, this validates the choice of the PML included in 
the FEM calculations. 

2. NUMERICAL METHODS 
2.1. Green’s dyadic technique 

We consider an object immersed in a homogeneous background medium. The addition of a sub¬ 
strate will be discussed later. The object and background are assumed nonmagnetic (/i = 1) with 
a relative permittivity, e and e^, respectively. In the following, we assume an exp(—icjt) time 
harmonic dependence for the fields. The Green’s dyad represents the electromagnetic response to 
an elementary excitation. Physically, the electric field scattered at the position r in presence of a 
point-like dipolar source po located at fq follows 

E(r) = ^G(r,ro)-po, (1) 

^0 

where cq is the free-space permittivity, ko the free-space wavenumber and G(r,ro) is the Green’s 
tensor associated with the whole system. It is an intrinsic electromagnetic quantity of the object, 
independent on the excitation process. For instance the local density of mode is given by 12 a ES] 

In addition, if the object is excited by an incident electric field Eq, the electric field can be expressed 
everywhere in the system by a 3D integral 

E(r) = Eo(r) + [[[ G(r, r') • AeEo(r')rfr', (3) 

J J Jobject 

with Ae = 6(r) — e^. Finally, the main numerical task is the evaluation of the tensor G for every 
couple of points (r, r'). The Green’s tensor can be computed by solving the self-consistent Dyson’s 
equation 

G(r,r') = Go(r,r')+ /c^ /// Go(r, r") • Ae • G(r", r')*" , (4) 

J J J object 

where Gq is the free-space Green’s dyad, that is analytical. In order to solve Dyson’s equation Q, 
the object is discretized into N cells of volume Vk (fc = 1, • • • , A^); 

N 

G(r,r') = Go(r,r') + fcp ^ Go"*(r, rk) • Ae • G(rk, r'), with (5) 

k=l 

Gjr*(r,rk)= flf Go(r,r")dr". (6) 

J J JVk 

Depending on the object shape, different meshes can be used. In case of spherical cells, the integral 
of the Green’s tensor over a sphere of radius centered at is analytical and writes m 


Go"*(r,rk) = C'fcGo(r,rk),with 

Airak ( sin{kBak) , A 

“ IT "‘“7 


, where ks — \f^kQ . 


(7) 

( 8 ) 
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Ck is a geometrical factor that reduces to the sphere volume Ck = 47ra|/3 = Vk of subwavelength 
size {kBdk 1). In this case, the discretized Dyson’s equationreduces to 13 HD] 

N 

G(r, r') = Go(r, r') + fcg ^ VfcGo(r, rk) • Ae • G(rk, r'). (9) 

k=i 

However, when plasmonic objects with high aspect ratio are considered, it is preferable to use 
elongated cells. And the volume integral Eq. has to be numerically evaluated. This is rather 
difficult because of the singularity of the Green’s tensor that occurs when the source and the 
observation points coincide (r' = r) j26j. In a smilar context, Chaumet et al used a Weyl expansion 
of the Green’s tensor to performed the numerical integration of Eq. in the reciprocical k-space 
m More recently, Massa and cowokers derived an approximate analytical expression for the 
polarisability of a cuboidal cell m- 

In the present work, we first transform the volume integral to a surface integral over the 
(cuboidal) cell surface before a numerical evaluation in direct space. This avoids the singular¬ 
ity at the cell center and convergence difficulties. 

In 2005, Gao et al demonstrated that this volume integral can be converted into a (flux) 
surface integral thanks to Ostrogradsky’s theorem m 


Go"*(r,rk) = [[[ Go(r,r")rfr" =-D(r) + // go(r, r")c?r". (10) 

J J JVk JJSk 

where D(r) = 1/A:^ if the observation point is within the integration volume (r G Vk) and is null 
elsewhere. For a rectangular parallelepiped mesh {ax b x c)^ centered at the point Fc = (xc, yc^ ^c)? 
the surface integral term writes for instance 


G 


int 

0,xx 


(r,rc) 


1 

47rfc^ 


ryc+bl2 

Jyc-b/2 


rzc-\-c/2 


{x — Xc — al2)e^^^^^^{ik\)Rx\ — 1 ) 




Gc-cI2 

{x - Xc + a/2)e'^^^^^^{iki,Rx2 - 1) 

^x2 


dzodyo . 


(11) 

( 12 ) 


with Rx 1^2 = [x — Xci: tt/2)^ -\- {y — yo)‘^ + ~ and ± refers to either Rxi (plus sign) or Rx 2 

(minus sign). Such surface integral is efficiently numerically computed using the Gauss-Kronrod 
method. The singularity of the Green’s dyad in the source region is properly taken into account 
in this approach. Another notable advantage of this approach is that it considerably reduces the 
memory cost for evaluating the Green’s dyadic. Finally, different mesh shapes can be considered by 
adapting the integral surface. And Dyson’s equation (Eq. is numerically solved using standard 
matrix inversion techniques. 

2.2. Finite element method 

The finite element method is a common numerical tool to study problems related to electromag¬ 
netism. Its main advantage is that it can be applied to irregular geometries. Since a detailed 
description can be found in Ref. [28], only a brief introduction will be given here. In classical 
electromadynamics, FEM is typically used to solve problems governed by Maxwell’s equations that 
can be formulated as a vector wave equation for the electric field 


V A [- V A E] - kocE = 0 . 

fi 


(13) 


The finite element method is based on the discretization of the whole geometry, including the 
surrounding medium, by simple elements (such as triangles in two dimensional problems, or tetra¬ 
hedral elements in three dimensional problems). The system of equations to be solved can be 
assembled after obtaining the weak formulation of the partial differential equations. 

However, one of the main difficulties in the finite element analysis of wave problems in open space 
is to truncate the unbounded domain. A common approach is to introduce a special layer of finite 
thickness surrounding the region of interest such that it is non-reflecting and completely absorbing 
for the waves entering this layer under any incidence. Such regions were introduced by Berenger 








4 


and are called perfectly matched layers (PML) [29] . The most natural way is to consider PML 
as a geometrical transformation that leads to equivalent e and fi (that are complex, anisotropic, 
and inhomogeneous even if the original ones were real, isotropic, and homogeneous) |30t. |3T] . This 
leads automatically to an equivalent medium with the same impedance than the one of the initial 
ambient medium since e and fi are transformed in the same way. The equivalent e and /i ensure 
that the interface with the layer is non-reflecting. 

3. RESULTS AND DISCUSSION 


L 



.P 


Figure 1: Benchmark problem: a silver nanorod of aspect ratio L/D in air. 

In the following, we consider a silver nanorod that consists of a cylinder rod capped with 
hemispherical ends (see Fig. [^. The diameter of the nanorod is D = 2R and the total length is L 
so that its aspect ratio is L/D. In this part of our numerical model, the silver nanorod is immersed 
in air. We use the silver dielectric function tabulated by Johnson and Christy [3^ . 



Figure 2: a) Near-held intensity spectrum calculated at one extremity (point P in Fig. located 4 nm away 
from the silver rod tip). FEM: Finite Element Method. AR = 5.33 and AR = 1 refer to GDT computation 
with cell aspect ratio a/b = 5.33 and a/b = 1 {b = 2.5 nm), respectively, b and c) Near-held intensity 
distribution maps calculated 5 nm above the rod surface at the two resonance peaks. Calculations were 
performed using hnite element method. The incident held is transverse magnetic (TM) and the incident 
angle is hxed at 7r/4. The intensity is normalized with respect to the incident intensity. The silver rod 
position is reported (white lines). 


We hrst investigate the near held optical response of the silver nanorod using FEM. The PML 
parameters (thickness and position) were determined by benchmarking the calculation of a well- 
known object characterized by the exact generalized Mie theory |33j (not shown). Considering 
dimer test structures, we achieve an excellent agreement with 50 nm thick spherical PML located 
at 150 nm from the nanostructure. We then consider the silver nanorod presented in Fig. A 
simple way to determine all the supported modes of the nanostructure is to excite it with an oblique 
incident plane wave Eq. We characterize the optical near-held response by computing the electric 
held intensity at one extremity (point P in Fig. [^. 
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Figure 2(a) [ displays the near-held intensity spectra of the silver nanorod (L = 100 nm and D = 
20 nm). We observe two surface plasmon polariton (SPP) resonances located around A = 440 nm 
and A = 682 nm. These resonances are in agreement with a Fabry-Perot resonator description 
[341 [35l[36]. The effective indices of the SPP guided along a 20 nm silver nanowire are neff = 2.76 
and neff = 3.98 at A = 682 nm and A = 440 nm, respectively. The cavity modes correspond to 
resonator lengths L -\r25 — mA/(2neff) with m = 1, 2,... m is the mode order and 6 refers to the 
field penetration depth into air. We obtain [m = 1, 5 = 12 nm] and [m = 2, 5 = 5 nm] for the two 
first mo des. Not e tha t these values for 6 reveal that high order modes are strongly confined [2]. 
Figures [2(b)| and 2(c) [ represent near-field intensity maps computed at these resonances located at 
A = 682 nm and A = 440 nm, respectively. The intensity map at A = 440 nm is not symmetric due 
to an interference with the incident excitation field (Fig. 2(c)[). The apparent symmetry observed 
at A = 682 nm is simp ly du e to the fact that the amplitude of the excited mode is stronger than 
the incident field (Fig. 2(b)[ ). 

We now apply the formalism of Green’s method presented in the previous section. The silver 
nanorod is discretized with rectangular parallelepipeds of dimensions a x b x c. We set b and c 
equal to 2.5 nm in the following, while varying the parameter a from a = 2.5 nm to a = 20 nm 
as shown in Fig. [^ The two hemispherical caps are discretized with a — b = c — 2.b nm for each 
case. Only the central cylindrical part of the nanorod is meshed with different a values. We use a 
transitory meshing area between the cap and the central part. The background is not discretized 
contrary to FEM calculations. 
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Figure 3: Near-field intensity distribution calculated 5 nm above the rod surface at the wavelength A = 
682 nm. The hemispherical caps are discretized with cubic cells of length 2.5 nm whereas the central part is 
discretized with rectangular parallelepipeds of dimensions axbx c with 6 = c = 2.5 nm and different aspect 
ratio AR — ajb. a varies from a = 2.5 nm (AR=1) to a = 20 nm (AR=8). Calculations performed using 
GDT. 


The near-field spectrum is calculated in Fig. [2(a)[ for two mesh aspect ratios, namely AR = 1 
(a = 2.5 nm) and AR — 5.33 (a = 13.3 nm). We observe an excellent agreement with the FEM 
calculated spectrum for the finer discretization, except for a peak near A = 730 nm. This peak 
disappears when changing the meshing size so that it is easily detected as a numerical artefact. 
In case of meshing aspect ratio AR = 5.33, the dipolar resonance (m = 1) is in agreement with 
the FEM computation. The m — 2 resonance is however slightly blue shifted. We represent the 
calculated intensity at A = 682 nm in Fig. [3 for mesh aspect ratio ranging from 1 to 8. When the 
aspect ratio of elementary cells is close to the one of the full object (Fig. [3(c)[ AR ^ 5) we can 
reach a very good agreement with the result obtained by the finite element method (Fig. 2(b)] ). 

Table indicates the number of meshing elements and the computing time for the different 
aspect ratios considered here. The computing time linearly increases with the cell number N 
since the limiting operation is the numerical integration of the free-space Green’s tensor over the 
cuboidal cell (instead of an increasing time as when the limiting factor is the resolution of the 
self consistent Dyson’s equation]^. 
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Aspect ratio 

1 

2 

5.33 

8 

Number of cells 

1944 

1112 

696 

488 

Computational time 

3hl3min 

lh43min 

lh04min 

48min 


Table 1: Summary of computational time and the number of cells discretization according to the aspect ratio 
of the elementary cell. Computations performed on a Intel X5650 (2.66 GHz) processor. 



Figure 4: Electric field intensity calculated 4 nm away from a single cell (blue curve, ’one cell’), a 3 x 3 x 3 cells 
object (red curve ’three cells’) and the complete structure (green curve). The cell aspect ratio is AR=5.33. 


Apart from the finest meshing, we o btain the better result for a meshing shape concomitant 
with the full object aspect ratio (Fig. |3(c)D . We attribute this to the fact that the resonance 
position of an elementary mesh is close to the object resonance. So the optical response of the 
whole structure is qualitatively described by the meshing and the object discretization refines the 
modelisation with a scaling law improvement. Figure represents the near-held intensity when 
considering a single cell (AR=5.33), 27 cells (AR=5.33) or the full object (AR=5). Although the 
localized plasmon resonance of a silver spherical or cubic nanoparticle is around A ~ 350 nm, the 
elongated cell presents a strong red shifted resonance around A = 920 nm and the resonance quickly 
converges to its hnal value as the object shape is built up. 

Finally, we characterize the held conhnement by computing the intensity as a function of distance 
to the rod’s extremity in hgure[^ FEM and GDT with high aspect ratio meshing are in quantitative 
agreement. The exponential ht shows a decay distance d = 11 nm in agreement with the value 
deduced from the Fabry-Perot resonator model ([m = 1 = 12 nm], see above). Note that the 

exponential decay is a simple way to characterize the mode conhnement but does not refer to the 
dipolar held decay [371 EH] • 


Figure [^represents the near-held intensity calculated at A = 440 nm for different meshing. We 
obtain a good agreement with the FEM-calculated map for the hnest meshing (AR=1 and AR=2). 
Larger mesh aspect ratios are too rough to grasp the fast spatial held variation over the rod length. 


The agreement between the FEM and the GDT m aps is not fnll> 
meshing (compare the maximum intensity in Figs. 


ips IS not iniiy qna n 
2(c)| and |6(a)|6(b)D . 


qua ntitative, even for the hnest 

_ _ Since the optical response 

is extremelly sensitive to the exact object shape, we attribute this diherence to the slight difference 
of the object shape due to the meshing. Indeed, GDT relies on a rectangular meshing so that 
we expect higher scattering on the edges. The corresponding object roughness is however of the 
order of rms ~ 1 nm so that the considered meshing is sufficient to describe object obtained within 
the state of the art of nanofabrication techniques. Considering a radius of 11 nm instead of 10 
nm, we obtain a quantitative agreement between the FEM (i? = 11 nm) and GDT (R = 10 nm) 
calculated intensities. We therefore conclude to the agreement between the two methods, within 
the nanofabrication tolerance. 
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Figure 5: Normalized near-field intensity calculated as a function of the distance to the rod extremity at 
A = 682 nm using FEM or GDT. For GDT, the hemispherical caps are discretized with cubic cell of length 
2.5 nm whereas the central part is discretized with rectangular parallelepipeds of dimensions a x b x c with 
6 = c = 2.5 nm and aspect ratio AR = a/b = 5.33. The green curve is an exponential fit I{d) = 
with A = 4653 and 5 = 11 nm. 
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Figure 6: Near-held intensity distribution calculated 5 nm above the rod, at the resonance of the second 
mode (A = 440nm). Galculations performed using GDT with meshes of different aspect ratios as indicated 
in the hgures. 


4. LDOS 


So far, we discussed the optical response of the metal nanorod to a plane wave excitation. This 
allows to detect the plasmon resonance and estimate the hied prohle of the mode. However, LDOS 
map better describes the mode since it is an intrinsic properties of the system. The LDOS is also a 
key quantity to interprete electron energy loss spectroscopy (EELS) [39l 00] and governs the decay 
rate of a quantum emitter in conhned system EH- It can also be manipulated at the nanoscale, 
opening promising perspectives to realize original nano-optical devices m- 

Figure [^represents the LDOS calculated above the nanorod using Eq. [^at the two resonances 
wavelength. As expected, the LDOS maps are symmetric and characterizes the suported modes, 
with a number m of node in agreement with the Eabry-Perot resonator description [l3]. The 
normalized LDOS amplitude of the hrst mode (up to 1700) i s com parable to the normalized 
electric held intensity calculated for a plane wave excitation (Eig. 3(a)| up to 1400). This again 
reveals the efficient excitation of this mode by a plane wave. On the contrary, the second mode is 
weakly excited by a plane wave, even at oblique incidence (compare the LDOS magnitude, up to 
PS 600 in Fig. and the electric intensity amplitudes, about ten times less in Figs. and[^). 
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Figure 7: Total LDOS map calculated 5 nm above the rod surface at the two resonance peaks. Calculations 
performed using GDT. LDOS is normalized with respect to its free-space value. 


5. EFFECT OF THE SUBSTRATE 

Finally, we discuss the effect of the substrate on the optical properties. This is an important 
point since plasmonic nanostructures are generally supported on a substrate. The presence of the 
substrate red shifts the resonance peak mi so that it has to be taken into account. However, it could 
lead to difficult FEM numerical implementation since the field radiatively leaks into the substrate 
(most refringent medium) so that it has to be sufficiently discretized and requests important memory 
resources. On the opposite, GDT easily takes into account the presence of the substrate without 
increasing the memory cost and with a minor increase of the computing time. In the present case 
where the Green’s dyad has to be integrated over the mesh surface, it is advantageous to work within 
the dipole image (quasi-static) approximation that leads to an excellent agreement with the exact 
retarded description for nanostructures laying on a glass substrate [H]. Let us consider a dipole 
Po = {Px^Py^Pz) at position ro = (xq, 2/o, ^o) above a substrate of dielectric constant Csub- The effect 
of the substrate on the dipole radiation is described by an image dipole at = {xQ^yo^ —^o) 
but in the homogeneous background of dielectric constant e^. If the observation point r is above 
the substrate, the image dipole to consider writes pim = {^sub~^B)/{^sub~\~^B){—Px^ ~Py^Pz)- If fhe 
observation point is below the substrate, the image dipole to consider writes = ‘2‘^sub/i^sub + 
^s)po- Since the Green’s dyad expresses the field scattered by a dipolar source, the integrated 
formulation is easily deduced within the image dipole approximation. For instance. 


/^int 
^ sub,XX 


(r,ro) = G'*o”‘^(r,ro) 


^sub ^B , 

^sub T ^B 


(^int ( 


(14) 


for an observation point r above the substrate. 

The near-held spectrum of the rod deposited on a glass substrate is calculated in hgure using 
GDT. We observe a strong shift of the resonances to the red. The resonance shift from A = 682 nm 
to A = 770 nm for the fundamental mode and A = 440 nm to A = 476 nm for the second order 
mode. The corresponding LDOS maps are represented in hgure We check that the mode prohles 
still correspond to the m=l and m=2 SPP modes. 

We also plot the near-held spectrum calculated using the FEM (Fig. |^. We observe good 
agreement with our GDT calculations except a lower intensity. We again attribute this to the 
extremelly sensitivity of the optical near-held response to the object shape. By considering a rod 
radius R=ll nm, FEM calculations lead to intensities comparable with GDT data. In addition, we 
would like to mention that the introduction of the PML layer is very critical for FEM in presence 
of the substrate. We observe strong variations of the calculated intensity with the PML center 
position. Best results are obtained when centering the PML at the rod (scatterer) center. If 
the PML is centered e.g. at the glass/air interface (that corresponds to a 10 nm shift only), the 











9 



Figure 8: Near-field intensity spectrum calculated at one silver rod extremity (point P, located 4 nm away 
from the silver rod tip). FEM: Finite Element Method. AR = 5.33 and AR = 1: GDT with cell aspect ratio 
a/b = 5.33 and a/h = 1 (with b = c = 1.6 nm), respectively. The rod is on a glass substrate of dielectric 
constant Csub = 2.25. The dotted curve represents the spectrum in absence of the substrate for comparison. 


calculated intensity is divided by 2 (not shown). Therefore, comparison with the GDT allows to 
determine optimized PML parameters (size and position) for FEM before considering more complex 
shapes. 
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Figure 9: Total LDOS map calculated 5 nm above the rod surface at the two resonance peaks. Calculations 
performed using GDT. The rod is on a glass substrate of dielectric constant Csub = 2.25. 


6. CONCLUSION 

In summary we used a dyadic Green technique based on a cuboidal meshing to compute the electric 
near-held and the LDOS near a silver nanorod. The fast variations of the electric held in the metallic 
nanostructures are well described by a numerical integration of the Green’s tensor over elongated 
meshes. Moreover, the computing efforts are reduced by transforming the volume integration over a 
mesh to a surface integral. The efficiency and accuracy of this method are verihed by comparing our 
results to those obtained by using hnite element method. This is a conhrmation of our technique, 
but reciprocically helps to determine the best parameters for FEM (notably PML). Green’s dyad 
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technique presented in this paper is very convenient, especially in the case of elongated structure 
(with high aspect ratio). Moreover, since the meshing has to be increased at positions where the 
field fastly varies, this method could be advantageously combined with top-down extended meshing 
algorithms [46] . In addition, the presence of a substrate is easily taken into account, demonstrating 
the versatility and efficiency of the method. Since the method relies on an numerical integration 
of the Green’s dyad over a mesh surface, more complex meshing (e.g. tetrahedral) could also be 
considered, allowing for a better description of the object shapes. Finally, as soon as the Green’s 
tensor of the whole structure is calculated, every electromagnetic responses of the system are easily 
determined (electric and magnetic near and far field, LDOS, optical forces, EELS response, ...). 
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